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Topological superconductors possess a nodeless superconducting gap in the bulk and gapless zero 
energy modes, known as “Majorana zero modes”, at the boundary of a finite system. In this work, 
we introduce a new class of topological superconductors, which are protected by nonsymmorphic 
crystalline symmetry and thus dubbed “topological nonsymmorphic crystalline superconductors”. 
We construct an explicit Bogoliubov-de Gennes type of model for this superconducting phase in the 
D class and show how Majorana zero modes in this model are protected by glide plane symmetry. 
Furthermore, we generalize the classification of topological nonsymmorphic crystalline supercon¬ 
ductors to the classes with time reversal symmetry, including the Dill and BDI classes, in two 
dimensions. Our theory provides a guidance to search for new topological superconducting materi¬ 
als with nonsymmorphic crystal structures. 
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I. INTRODUCTION 

The research on topological superconductors (TSCs) 
has attracted intensive interests due to its gapless 
boundary excitations, known as the “Majorana zero 
modes’’^"— , with intrinsically non-local nature and ex¬ 
otic exchange statistics, and aims in the potential ap¬ 
plications in low-decoherence quantum information pro¬ 
cessing and topological quantum computation&i^^— . The 
search for new topological superconducting phases and 
materials is a substantial step for this goal. 

The first classification of TSCs (and also other topo¬ 
logical insulating phases) was achieved by Schnyder et 
al^ based on Altland-Zirnbauer symmetry class^i^ 
for the systems with or without particle-hole symmetry 
(PHS), time reversal symmetry(TRS) and their combi¬ 
nation, the so-called chiral symmetry. Later, it was real¬ 
ized that when additional symmetry exists in a system, 
new topological phases can be obtained, and the gap¬ 
less edge/surface modes require the protection from ad¬ 
ditional symmetry. In particular, it has been shown that 
new topological insulating and superconducting phases 
emerge when the system has mirror symmetry and are 
dubbed “topological mirror insulators’— and “topo¬ 
logical mirror superconductors” ^ respectively. Re¬ 

cent work has also revealed that nonsymmorphic sym¬ 
metry, including glide plane symmetry and screw axis 
symmetry, can lead to new topological insulating phases, 
as well as topological semi-metal phases^^—. In this 
work, we are interested in the role of the nonsymmor¬ 
phic crystalline symmetry, mainly glide plane symme¬ 
try, in the classification of TSCs. We focus on the fol¬ 
lowing three questions: (1) are there any topologically 
non-trivial phases that are protected by glide plane sym¬ 
metry? (2) What’s the difference between glide plane 
symmetry and mirror symmetry in the classification of 
TSCs? (3) What’s the relationship between this super¬ 
conducting phase and other TSCs? Below, we will first 
discuss the role of glide plane symmetry in the classifi¬ 


cation of superconducting gap functions, which indicates 
the possibility of topological superconductors protected 
by glide plane symmetry, thus dubbed “topological non¬ 
symmorphic crystalline superconductors (TNCSc)”. We 
also construct an explicit tight-binding model in the D 
class with boundary Majorana zero modes and demon¬ 
strate that the existence of Majorana zero modes comes 
from nonsymmorphic symmetry of this model. Finally, 
we discuss the relationship between TNSCs and weak 
TSCs and generalize TNCSc to the classes Dill and BDI 
with time reversal symmetry for both spinless and spin-^ 
fermions. 


II. NONSYMMORPHIC SYMMETRY AND 
SUPERCONDUCTING GAP FUNCTION 

In this section, we will first consider the role of glide 
plane symmetry in the classification of superconducting 
gap functions. We start from a generic Bogoliubov-de 
Gennes (BdG) type of Hamiltonian of superconductors 
with nonsymmorphic symmetry in the normal states, 
which can be written in the momentum space as 

^ ^ c\)HBdG ^ 


with 


HbcIG 


h(k) - A(k) \ 
At(k) -h*{-k) + fij 


( 1 ) 


where h(k) is for single-particle Hamiltonian of nor¬ 
mal states, /i is the chemical potential and A denotes 
the superconducting gap function. Ck is an annihila¬ 
tion operator with n components and we also use Ck,a 
(a = I,...,n) to denote each component with a = {s,/} 
for spins s and orbitals (lattice sites) 1. The supercon¬ 
ducting gap function is related to annihilation operators 
by AQ,^/ 3 (k) = Vb{ck,/ 3 C-k,a), where Vb is the strength of 
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attractive interactions. The BdG Hamiltonian satisfies 
the PHS CHBdG(}^)C~^ = —HBdG{—^) with the PHS 
operator (7 = ri x ITsT, where ti is the first Pauli matrix 
acting on the Nambu space, I is an n x n unit matrix and 
K is complex conjugation. The PHS (or Fermi statistics) 
requires the constraint A(k) = —A^(—k) for the gap 
function. 

Next, we consider how the nonsymmorphic symmetry 
yields constraint on the forms of single-particle Hamil¬ 
tonian and superconducting gap functions in a nonsym¬ 
morphic crystal. Here we consider the glide plane sym¬ 
metry, represented by g = {m|r} where m is a mirror 
operator and r is a non-primitive translation operator 
along a direction within the mirror plane. For single¬ 
particle Hamiltonian, the glide plane symmetry requires 
£)^(g)/i(k)£)k( 5 ) = where D\^{g) is the representa¬ 

tion matrix for glide plane symmetry at the momentum k 
and defined as = E /3 ^k.a/ 3 ( 5 )Cgk./?^' ^ere we 

emphasize that the representation matrix for glide plane 
symmetry takes the form Dk{g) = D{m), where 

is a phase factor due to a non-primitive translation and 
D{m) is the projective representation of mirror operator 
TO. For the case with only glide plane symmetry, all the 
projective representations are one dimensional (ID) and 
equivalent to the conventional representations. 

The symmetry of the gap function A(k) is deter¬ 
mined by the Cooper pair wave functions, which trans¬ 
form as the direct product of the representation Dl^{g) (g) 
^lk(5)^- For the case with only glide plane symmetry, 
all the ID representations can be labeled by Dk)^) = 
e*k r£)( 77 jj _ ijgjk-r ^ _ j-j fgj. gpin-i systems 

and J = ±1 for spinless systems. Thus, the gap func¬ 
tion should transform as Dl.{g)A{k)D'^y^{g) = T]A{gk), 
where 77 = ± applies for both the spin-i and spinless sys¬ 
tems and depends on the nature of superconducting gap 
functionals. We will show how to classify different super¬ 
conducting gap functions based on glide plane symmetry 
explicitly for a model Hamiltonian in the next section. 

We emphasize that the superconducting gap func¬ 
tion may preserve (ry = -f) or spontaneously break 
(?7 = —) glide plane symmetry. Nevertheless, sim¬ 
ilar to the case of inversion symmetrySiiSS or mirror 
symmetry^, one can always re-define a glide plane sym¬ 
metry operation as G,,(k) = Diag[Di^{g),r]D'^^^{g)] for 
the BdG type of Hamiltonian, which satisfies the condi¬ 
tion G~^{k)HBdG{^)Gr,{k) = HBdGig^)- In this way, 
we can regard the BdG Hamiltonian as a semiconductor 
Hamiltonian with additional PHS. 

Due to the existence of the glide plane symmetry 
G,,(k), the eigenstates '^(k) of the BdG Hamiltonian, 
HBdG'^{^) = Etl^ik), can also be chosen to be the 
eigenstate of G^(k), Grj{k)ip{k) = ^^e*’^'^ 7 /^(k), on the 
glide invariant plane (GIP) in the momentum space, 
(/k = k (mod P), where P is a reciprocal lattice 
vector. Here <5^ is given by ± (±j) for the spin¬ 
less (spin-i) systems and we call the eigenvalue 
as glide parity. Next, we look at the relationship 
of glide parities between one eigenstate ^/’(k) and its 


(a) G+ configuration (b) G- configuration 



FIG. 1. (Color online). Two different configurations for 
G±(k). For G+(G_), the Hamiltonian symmetry class is 
D(A) along k ■ r = 0 lines; while along k • r = -I lines, the 
Hamiltonian symmetry class is A(D). The red dots denote 
possible Majarona zero modes at ends of the lines. 

partner '!/!’(—k) = Cip{k) under PHS. Direct calcula¬ 
tion gives G^(—k)7/^(—k) = r]S*e~^^''^tl;{—k) by using 
that GG,,(k)G“^ = r]Gn{—k)^. Therefore, i/^k and its 
particle-hole partner possess glide parity and 

77(5*6“®*''®’, respectively. This leads to the conclusion as 
depicted in Fig. [T] When the gap function satisfies 
G-i-(k) symmetry, for the spinless (spin-i) systems, ipk 
and its particle-hole partner tp-k share the same glide 
parity along the momentum line k-T = 0(k-r=^)on 
the GIP, while they have opposite glide parities along the 
momentum line k • r = ^ (k • r = 0) on the GIP. When 
the gap function satisfies G_(k) symmetry, we find an 
opposite behavior for the momentum lines k • r = 0 and 
k • r = I, compared to the case of G+ (k) symmetry. We 
notice that the momentum line k ■ r = ^ corresponds to 
the BZ boundary since 2t is a primitive lattice vector of 
the system. 

Here we emphasize different roles of glide plane sym¬ 
metry and mirror symmetry for the BdG Hamiltonian 
of superconductivity. For the glide plane symmetry 
g = {toIt} and the corresponding mirror symmetry to, 
the GIP and the mirror invariant plane are the same. As 
shown in Ref. [ 2 ^, the PHS either preserves the subspace 
with a fixed mirror parity or transforms the subspace 
with one mirror parity to the other. In contrast, due to 
the additional phase factor from the non-primitive trans¬ 
lation of glide plane symmetry, the behaviors of PHS act¬ 
ing on the glide parity subspaces are always opposite for 
the momentum lines k-r = 0 and k-r = This prevents 
us to define a topological invariant on the whole 2D GIP 
since two glide parity subspaces are always “connected” 
to each other. However, if we limit the glide parity sub¬ 
space only on the momentum line k-r = 0ork-T = ^, 
the PHS will either preserve the glide parity subspace 
or transform the subspace with one glide parity to the 
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other, similar to the case of mirror symmetry. This im¬ 
mediately suggests the possibility of defining topological 
invariants on the ID momentum lines k-r = 0 or k-r = ^ 
for superconductors with glide plane symmetry. Below, 
we will present explicitly a BdG type of model Hamilto¬ 
nian with glide plane symmetry and show the existence 
of Majorana zero modes at the boundary. Then we will 
discuss bulk topological invariants and the corresponding 
topological classification. 


III. MODEL HAMILTONIAN IN THE D CLASS 

Our spinless fermion model with glide plane symmetry 
is based on a two dimensional (2D) rectangle lattice with 
two sets of equivalent sites, as shown by A and B sites 
in Fig. UKa) and (b). The glide plane symmetry operator 
is given by gz = {m^lr = (1,0,0)} with a reflection 
along the z direction followed by a translation of a/2 
along the x direction (a is a lattice constant), and relates 
the A sites to the B sites. The normal state Hamiltonian 
reads 

h{k) =e(k)(To -f t 3 Cos(^^^^Y^)cos(^)cri 

, Akx ~'P)^\ ■ /k^a 

+t3Cos( - - - )stn[—)a2 ( 2 ) 

on the basis |A,k) and \B,k), where e(k) = mo + 
ticos{kxa) + t 2 Cos{kya), (Tq is a 2x2 unit matrix, ai 
with i = 1, 2, 3 are Pauli matrices that describe 
the A and B sites and (j) depends on the choice of 
orbitals^. Furthermore, the glide plane symmetry op- 
erator on such a basis is Dk{g) = (cos(^)(Ti -b 
sin{^^)(T 2 ). One can easily check that D^(g) = 
and D^\g)H{kx,ky)DAg) = H{kx,ky). 

As discussed above, the gap functions can be classified 
according to glide plane symmetry and when the glide 
plane symmetry for the BdG Hamiltonian is G^, the gap 
function satisfies three conditions: A^(k) = — A(—k) 
(PHS); Dl^{g)I\(k)D*_^{g) = g/2^{gk) (glide plane sym¬ 
metry) and A(k) = A(k + G)^. The complete classi- 
ficaiton of gap functions for this model Hamiltonian is 
discussed in the Supplemental Material^. Here we only 
consider two typical gap functions A+ = Aosin{kya)ao 
and A_ = Aosin{kya)a 3 with the symmetries G+ and 
G_, respectively. We take the BdG Hamiltonian (Eq. 
[T|) with the single-particle Hamiltonian (Eq. [2|) and the 
gap function A± and calculate energy dispersion of this 
Hamiltonian on a slab configuration. The slab is chosen 
to be infinite along the x direction and finite along the 
y direction, so that the glide plane symmetry gz is still 
preserved. The energy dispersion is shown in Fig. [2] (c) 
for A+ and (d) for A_. In both cases, one can find two 
edge bands appearing in the bulk superconducting gap at 
one edge. However, these two edge bands cross at zero 
energy and give rise to Majorana zero modes at F for A+ 
(Fig. HKc)), but at A for A_ (Fig. [^d)). 



FIG. 2. (Color online), (a) and (b). Schematic plots of the 
lattice structure from top view and side view. They are ID 
chains along y direction. There are two inequivalent atom 
sites, denoted as A(Red ball) and B(Black ball), respectively. 
A plane passing through the dashed green line is the glide 
plane, (c) Edge modes for G+ configuration with A+. (d) 
Edge modes for G- configuration with A_. (e) Brillouin 

zone(Black square) and extended Brillouin zone(Red dashed 
rectangle) defined by glide plane symmetry, (f) A general 
dispersion for a ID chain with glide plane symmetry. 


The underlying physical reason of different positions 
of Majorana zero modes for these two cases comes from 
the relation between glide plane symmetry and PHS dis¬ 
cussed in the last section. Let’s take the case of the 
gap function A_|_ with G+ symmetry as an example. The 
state '0k and its particle-hole partner ip-k. share the same 
glide parity at F {k^ = 0), and thus it is possible for them 
to be the same state. Since PHS changes the energy E 
of '0k to —E of 0-k, the eigen energy must be zero once 
they are the same state. This analysis also suggests that 
two Majorana zero modes at F must belong to different 
glide parity subspace, and thus no coupling is allowed 
between them to open a gap. In contrast, the glide par¬ 
ities for 0k and tp-k are opposite at X {kx = f )■ Thus, 
these two states must be different at X and PHS can not 
require their energies to be zero. This analysis can also 
be applied to A_ with G_ symmetry and leads to the 
opposite conclusion. Another intuitive picture to prove 
non-trivial properties of ID edge modes in Fig. [21(c) and 
(d) is to consider a general one dimensional superconduc¬ 
tor with glide plane symmetry. As shown in Ref. 
due to the glide plane symmetry, all the bands must ap¬ 
pear in pairs, as shown schematically by two black lines 
(two bands with opposite glide parities) in Fig. [2] (f). 
Furthermore, the PHS of superconductivity requires two 
additional hole bands at the negative energy, as shown 
by two red lines in Fig. |2](f). Therefore, there must be 
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even number of pairs of bands for a ID nonsymmorphic 
superconductor. A single pair of bands shown in Fig. [2] 
(c) and (d) can only exist at the ID boundary of a 2D sys¬ 
tem. This gives the “no-go” theorem for nonsymmorphic 
superconductors^. 


IV. BULK TOPOLOGICAL INVARIANTS AND 
THE EXTENDED BRILLOUIN ZONE 

The above analysis has shown that two momentum 
lines k • r = 0 and k • r = ^ play the essential role in the 
classification of TSCs in nonsymmorphic crystals. We 
can view the bulk Hamiltonian onk-T = 0ork-T = ^ 
as a ID Hamiltonian. For the case of A+ with the G+ 
symmetry, the Hamiltonian HBdG (Eq- Hj) has PHS along 
the line k • r = 0 for each glide parity subspace, thus be¬ 
longing to the D class, while it has no PHS along the line 
k • T = ^ for each glide parity subspace, as shown in Fig. 
□ (a). Since two glide parity subspaces are decoupled, one 
Z 2 topological invariant of the D class can be defined on 
the line k • r = 0 in the glide parity subspace for a ID 
Hamiltonian. In contrast, for the case of A_ with the G_ 
symmetry, one Z 2 topological invariant can be defined on 
the line k • r = ^. In our example, we can re-write the 
BdG Hamiltonian with the eigenstates of G± as a basis 
and one can see immediately for the case with the G+ 
(G_) symmetry, the Hamiltonian is exactly equivalent 
to the ID Kitaev model of p-wave superconductors^ in 
each glide parity subspace when kx = 0 {kx = f }^- 

More insights about this system can be obtained from 
the view of the extended Brillouin zone (BZ)^, which has 
been widely used in the field of iron pnictide supercon¬ 
ductors. For nonsymmorphic crystals, all the eigenstates 
of the Hamiltonian can be labeled by the eigenvalues of 
glide operators, defined as G^'0(k) = e®'^''”^(k), in which 
k is called “pseudocrystal momentum”— and defines the 
extended BZ. For our model, the glide plane symmetry 
operation only involves translation by | along the x di¬ 
rection, and thus the extended BZ for kx is doubled along 
the X direction {kx & compared to the con¬ 

ventional BZ for kx, as shown in Fig. [2] (e). Since we 
have G,,'!/'(E) = Srje^^''^'tlj{'k), this suggests that the pseu¬ 
docrystal momentum k is related to momentum k by 
k = k when Sr/ = + and k = k-l-Q with Q • r = ±7r when 
Srj = —. Here the sign of Q • r is determined by keeping 
kx in the region and ky in the region f ]■ 

As a result, the BdG Hamiltonian can also be rewrit¬ 
ten as i7f“^g(k) = HBdG,+ {^) = HBdG,+ {^) for kx S 
[-f, f] and H%-^^{k) = HBdG,-{^) = i^BdG.-(k - Q) 
for kx G [-,—1 and kx & — in the extended 

BZ. Here HBdG,± is the BdG Hamiltonian in the sub¬ 
space with glide parity For our model Hamil¬ 

tonian, HBdG,± corresponds to the two by two Hamil¬ 
tonian defined in the Supplemental Material^^. For the 
case of A+, the form of is given by iJg^Q(k) = 

— {e{k)—fd+t3COs{^)cos{^))T3+t3sin{^)sin{^)To+ 


AQsin{kya)Ti, where e(k) = mo+ticos{kxa)+t 2 cos{kya). 
We notice that if we take the hopping parameters ti and 
ts along the x direction to be zero, this Hamiltonian ex¬ 
actly corresponds to the ID Kitaev chain with one Majo- 
rana zero mode at the open boundary^. With the hop¬ 
ping along the x direction, all the ID Kitaev chains are 
coupled along the x direction, so Majorana zero modes at 
the end of the chains couple to each other and expand into 
a band. This corresponds to the weak TSGs^i^, which 
is in analogy to weak topological insulators^. The PHS 
requires A(k) = — A(—k) for the band of Majorana zero 
modes. Therefore, zero energy states can only appear for 
kx = 0 and kx = — {kx is periodic in —), which both 
correspond io kx = Q in the conventional BZ. In contrast, 
for the case of A_, the gap function comes from the so- 
called 77 pairing for two electrons with the momenta k and 
Q — k to form a Cooper pai r^^d^- ^a (Q = (^,0) for our 
model). In this case, the PHS requires A(k) = —E{Q—k) 
for the Majorana band, leading to the zero energy states 
a,t kx = This analysis based on the extended BZ is 
consistent with our previous results and show explicitly 
the relationship between TNSCs and weak TSCs. 


V. DISCUSSION AND CONCLUSION 


The above results for TNCSc can be directly gener¬ 
alized to the systems with spin-i and with additional 
time reversal (TR) symmetry. For spin-i systems, since 
5y in the glide parity is given by ±i, there is an addi¬ 
tional minus sign when considering how the glide parity 
of an eigenstate of the BdG Hamiltonian transforms un¬ 
der PHS. This leads to the consequence that the Z 2 topo¬ 
logical invariant can be defined at k • r = ^ (k • r = 0) 
for the systems with the G-). (G_) symmetry. According 
to the standard topological classification, TR symmetry 
can change the symmetry class from the D class to BDI 
for spinless systems and Dill for spin-i systems. To see 
how it affects the classification of TNCSc, we consider an 
example of a spin-i system in the Dill class with the 
G_|_ symmetry. If we take a state ^/’(k) with glide parity 
where 5 = i, the glide parity of its PHS partner has 
been shown to be ” and the glide parity of its TR 

partner is also , where 0G+(k)0“^ = G+(—k) 


is used and TR operator is 0 



with 


0e = iaQS 2 K and S 2 the second Pauli matrix acting on 
spin space. One can see that chiral symmetry H = G x 0 
exists in each glide parity subspace for any momentum. 
In addition, at the momentum line k ■ r = PHS and 
TRS also exist in each glide parity subspace. There¬ 
fore, the symmetry class is Dili for the momentum line 
k ■ r = ^ and AHI for other momentum lines (k • r 7 ^ 
in each glide parity subspace. This leads to Z 20 Z 2 clas¬ 
sification at k • T = ^ , Z classification at k • t = 0 and 
Z 0 Z classification at other momentum lines for the whole 
BdG Hamiltonian2£i^, in sharp contrast to the Z x Z 
classification of topological mirror superconductors in the 
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Dill class^^. This classification leads to the existence of 
edge flat bands in the Dill class for TNCSc (See Supple¬ 
mental materials^). The spinless and spin-i TNCSc in 
classes D, Dill and BDI are also studied in the Supple¬ 
mental Material2£. Nonsymmorphic symmetry is known 
to exist in several classes of superconducting materi¬ 
als, including iron pnictide superconductors^^— , BiS 2 - 
based layered superconductors^^— , and heavy fermion 
superconductors^ii^, e.g. UPta^, UBeig^. Our topo¬ 
logical classification of TNCSc can be directly applied to 
these systems to search for realistic topological supercon¬ 
ducting materials. 
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Note added. - After finishing this paper, we notice 
a paper on arxiv^, which concerns possible topological 
superconducting phases in monolayer FeSe and poten¬ 
tial relation to nonsymmorphic symmetry. We also no¬ 
tice another recent paper on arxiv2S about topological 
classification of TNCSc based on the twisted equivariant 
K-theory. 


Appendix A: Bogoliubov-de Gennes Hamiltonian and the glide plane symmetry 


We start from a generic mean-field Hamiltonian of superconductors with nonsymmorphic symmetry in the normal 
states, which can be written in the momentum space as 

~ M)Ck,/3 + ■^(^a,/3''-k.aCk,/3 + 


k,CK,/3 


k,a,/3 


= ( JT 


"-k 


(Al) 


with 


HBdG 


h(k) - A(k) \ 
At(k) -h*{-k)+^lJ 


(A2) 


where h(k) is for single-particle Hamiltonian of normal states, /r is the chemical potential and A denotes the super¬ 
conducting gap function. Ck is an annihilation operator with n components and we also use Ck,Q; (ci = l,...,n) to 
denote each component with a = {s,/} for spins s and orbitals (lattice sites) 1. The superconducting gap function is 
related to annihilation operators by AQ,,g(k) = Vo(ck,,gc_k.a), where Vq is the strength of attractive interactions. The 
Bogoliubov-de gennes (BdG) Hamiltonian satisfies the particle-hole symmetry (PHS) CHBdG{k)C~^ = —HBdG{—k) 
with the PHS operator C = ri x lAT, where ti is Pauli matrix acting on the Nambu space, I is an n x n unit matrix 
and K is complex conjugation. The PHS (or Fermi statistics) requires the constraint A(k) = —A^(—k) for the gap 
function. 

glide plane symmetry can be expressed as g = {m\T} with m a mirror operator and r a non-primitive transla¬ 
tion operator along the direction within the mirror plane. We have api9)‘^l]i g ^^id gck_,ag~^ = 

p Dk,ap{g)cgk,p where D\^{g) is representation of space group g in the little group of wave-vector k. Further, 
Dk({m|T}) = e**^ '^D(m), where Dm is the projective representation of mirror operator m. 

Since the normal state a^«/ 3 (^^)ck ,/3 has the glide plane symmetry, gHeg~^ = He- Explicitly, 


^hap(k)gCk,ag ^k.a7(5)4k,7^a/3(k)^k,/3A(ff)Cgk.A 

k,CK/3 k,a,/3 

= 4,7(X1^9-'k.a7(ff)^a/3(5"^k)Dg-lk_^A(ff))ck,A 

k. 7 A a/3 


Therefore, we have 


^*g-^k,ai(9)Kp{g ^k)Dg-ik,,3A(ff) = h^A^) 


(A3) 
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More elegantly, the normal state Hamiltonian under glide plane symmetry satisfies 

Dl{g)h{k)D^{g) = h{gk). (A4) 

We also have Dl_{g) = D^^{g). 

The pairing terms in Eq. lAll fulfill the requirement of glide plane symmetry. However, the gap functions are not 
necessary to respect the glide plane symmetry. 

giY^ ^c.piy)cYcYp)9~^ = Y ^k.a7(ff)^-k./3A('?)Aa/3(k)cJk.aC-gk./3 = Z! ZA(k)cZclk,A 

k.a/3 k,a/37A k.7A 


where A.YA(k) is the transformed gap function. 


Thus, we arrive at A.yA(k) = 


Sa/3 ^g-ik a7(5)'^-g~ik ^k). This indicates that the gap function A(k) transforms according to 

the decomposition of the direct product of the representation Dl^{g) (g) D'^^{g). 

If we choose the basis as the eigenstates of glide plane symmetry operator, all representations are reduced to one 
dimension(lD) and Hklff) = D{m) = where 5 = for spin-^ systems and <5 = ±1 for spinless systems. 


Thus, A.yA iy)=D 


S'^'-k.77 


{9)D: 


5''^k,AA 


^k) = e 


:>-W 




‘k-r 


SlAapig ^k) = S*SlAap{g ^k). We can 


write down the requirement of gap function under nonsymmorphic symmetry in a compact way 

A( 5 k) == Dl{g)A{k)Dl^{g) = gAigk) 


(A5) 


where ry = ±1 applies for both spin-i and spinless systems. 

The matrix form of the BdG Hamiltonian is shown in Eq. IA2I in the basis of Nambu space dt(k) = 

According to Eq. IA51 we can treat the BdG Hamiltonian as a semiconductor Hamiltonian with additional PHS and 
re-define the glide plane symmetry operator as 



G,(k) 


Di.{g) 0 \ 

0 vD\{9) J 


(A6) 


Next we check how the BdG Hamiltonian transforms under glide plane symmetry by considering possible gap 
functions required in Eq. IA51 For the case with Dl.(g)A(k)DZi^(g) = gA(gk), we have 


G--(k)flB.a(k)G,(k)= ('‘aw" 

^ f H9k)-9 ^k^(ff)^(k)?7^1k(ff) 
h.c. —h*(—gk) + /i 

^ f Hgk) - 9 ’7^A(gk) \ 

h.c. —h*{—gk)+fi J 

^ ( K9^) - 9 A(gk) \ 

h.c. —h*{—gk)+g, j 

= HBdoig^) 


A(k) 

-h*(-k) + /i 


DUg) 0 \ 

0 vD*_^{9) j 


(A7) 


where Dl^{g) = D^^{g) is used. This gives us the form of symmetry transformation for the BdG Hamiltonian. 

Due to the glide plane symmetry, all the eigenstates of the BdG Hamiltonian at the glide invariant plane (GIP) 
can be also expressed as the eigenstates of glide plane symmetry and the corresponding eigenvalues are dubbed 
“glide parity”, as discussed in the main text. To show the glide parities of a state and its particle-hole partner, we 
take an example of the case with the G+ symmetry. If the gap function preserves the glide plane symmetry, the 
BdG Hamiltonian commutes with G+(k) on GIPs. This indicates that one can simultaneously block diagonalize 
HBdG and G+(k) with a set of common eigenvectors. Each block owns a glide parity with (5+ = ±l(±i) for 

spinless (spin-i) systems. 


We start from that GG+(k) = 


K 


( Dkig) 

V 0 


D 


0 

lk(5) 


D-u{g) 0 
0 D*^i9) 


0 1 
1 0 


K = G+(—k), i.e. 


GG+(k)G“^ = G+(—k), where G is the PHS operator. One can pick up a common eigenstate V^(k) oiHBdG and G+(k) 
with glide parity (5+e*''''”. Its particle-hole partner is denoted as V'(~k) = G^(k). Then we have G+(—k)^/;(—k) = 
G+(—k)G^(k) = G(G+(k)'!/;(k)) = G(5+e®*^'”'!/'(k)) = Therefore, ipk and its PHS partner ^_k 


possess glide parity and S^e 


respectively. 
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Appendix B: Model Hamiltonian 

In this section, we will show how to construct a tight-binding model for the TNCSc. We use (j){Y — R,; — r^) 
to denote Ldwding orbital a at r, where is the position of unit cell and is the position of atom site 
or orbital inside a unit cell. The Bloch wave function is dehned as ^/’k,a(r) = — R^ — r^) by 

using the linear combination of atomic orbitals(LCAO). It should be emphasized that the phase factor in our 

construction does not include the position Ya. This Bloch wave function can be written as '!/'k,c((r) = e*’' '’Mk,c((r) 
where Uk,a(r) = ^ It is easily checked that Uk,a(r + R) = X^r^ + 

R-Ri-r^) = -^ X^Ri -!■«) = + = Mk,a(r)- Thus, the 

Bloch theorem holds for such a choice of LCAO. Under such a choice of LCAO, we have that V'k,a(r) = '0k-i-p,c((r), 
where P is a reciprocal lattice vector. Any Hamiltonian on such a basis iJa,/ 3 (k) = / (ir'0^ j,^(r)iJ'0k,,s(r) satisfies that 
i 7 „,; 3 (k-pP) = /drV'k+p,c.(j')^^k-KP,/ 3 (i-) = /drV'k,„(r)HV’k,/ 3 (r) = Ha,p{k), i.e. 

i7(k-pP) = i7(k). (Bl) 

Similarly, for the gap function, A(k 4 . P) = (/ (irb'k+p,„(r)-0k-rp,a(r)) = (/ drV'k.a(r)-0k.a(i')) = A(k), i.e. 

A(k-t-P) = A(k) (B2) 


Next let us take cj^ ^ and Ck,ct(cJ ^ and Cr,a) to present creation and annihilation operators of 'i/;k,a(r)(<^(r — Ri — ra)) 
and show the real space form of the tight-binding model discussed in the section “Model Hamiltonian” of the main 
text. The normal state tight-binding Hamiltonian for the lattice structure reads 


Be =mo 






[ E 

i,s—A.B 


t 


t2 t 

2 i-\-dy 


[ y (4 .aC*.b + cl^Ci-d^^B) + B.C.] (B3) 


where i = denotes index of unit cells, dx,y are primitive lattice vectors along x and y direction, A and B 

denote two inequivalent atom sites and H.c. represents their conjugation parts. We further obtain a tight-binding 
model in the momentum space by performing an unusual Fourier transformation cj^ ^ s = 

e'^^'^Tk.s, where R* is the position of the i*^ unit cell. Such a Fourier transformation simplifies the Hamiltonian 
and leads to He(k -|- P) = He(k) with P is a reciprocal lattice vector. 

The gap functions for G± configurations mentioned previously need to satisfy three conditions: (1) PHS 

A^(k) = -A(-k); (B4) 

( 2 ) glide plane symmetry 

Di{g)A(k)Dl^,ig) = r;A(gk); (B5) 

and (3) 

A(k) = A(k-f G). (B 6 ) 

These three conditions allow us to classify all the possible gap function for this model Hamiltonian. Let us define 
the gap functions Ai(k) = /i(k)ri(k), where i = I,...,4, /i(k) is a complex function of k and ri(k) are four 2x2 
matrices, as shown in the second column in Table HI 

The second and third column in Table U shows the “parity” of F matrices in the sense of PHS and glide plane 
symmetry. Due to the PHS (IB4I) . the parity of /i(k) is determined by the parity of F from the second column, which 
is listed in the fourth column. For 2D system, k = {kx, ky) if g = {m^jr}. For the last two columns, the existence of 
fi{k) is determined by Eq. (IB5|) and in the third column. Let us take the r|i matrix with = I as an example. 
For G+, since glide operation does not act on /i(fc), we obtain /i(/c) = fi{k) and there is no constraint on fi. But 
for G_, Eq. (IB5I1 and = 1 together requires fi{k) = —/i(fc), leading to fi{k) = 0. Thus, no term is possible for 
G_. Similar analysis can be applied to other matrices. Based on the parity of /i(k) on the fourth column, we can get 
possible polynomials of /i(k), as listed in Table HIl 

The parameters for the calculation of energy dispersion of the BdG Hamiltonian are shown in Table IHII 


TABLE I. Here Cc, Cs and ^ are defined by rf(k) = ^cTii-k), (g)ri(k)DLk(5) = ? 9 ri(k), /i(k) = C//i(-k). PHS 
and glide plane symmetry requirements on fi(k) provide the = ±1 and = ±1 on the second and third columns. The 
= ±1 in the fourth column is the parity of /i(k) and obtained from the PHS (IB4II and ^c- NA in the fifth and sixth columns 
represents “not available”. Here Pi = cro, r 2 = ua, Pa = cos{^^)ai + sin(J^^)a 2 and r 4 = sin{j^)G\ — cos(^f^)cr 2 . 



^c| 


C/ 1 

G+: /i(k) 1 

G-: /,(k) 

Pi 

+ 

+ 

- 

Valid 

NA 

T 2 

+ 

- 

- 

NA 

Valid 

Ps 

+ 

+ 

- 

Valid 

NA 

P 4 

- 

- 

+ 

NA 

Valid 


TABLE 11. possible polynomials of /i(k) for each Ai(k). NAs in the fourth and fifth columns represent “not available” 




G+: /i(k) 

G-: /Gk) 

Ai(k) 

= /i(k)Pi 

sin{kxa), sin{kya) 

NA 

A2(k) 

= /2(k)P2 

NA 

sin(kxa), sin{kya) 

A3(k) 

= /3(k)P3 

sm(¥) 

NA 

A4(k) 

= /4(k)P4 

NA 

cos(%^) 


Appendix C: Hamiltonian in the extended Brillouin zone for G± 


In this section, we will analyze our model Hamiltonian in the extended Brillouin zone. For the case with the G+ 
symmetry, G+(k) = e*^(cos(^)ro 0 cti + sin{^^)To 0 0 - 2 ). The eigenvalues are =Fe*^- The eigenvectors are 
= ( 0 , 0 ,l)"^, U 2 - = (— 1 , 0 , 0 )^, = ( 0 , 0 ,e“*^, 1 )^ and U 2 ,+ = , 1 , 0 , 0 )"^, where 

Uzp corresponds to eigenvalues On the above eigenvectors, the BdG Hamiltonian in the case with the G+ 

symmetry can be expressed in a block diagonal matrix, which reads 




I — ,g+ 

I Aosin(kya) 
0 
0 


HBdG,g+ = 


Aosin{kya) 

y-,g+ 

0 

0 


0 0 

0 0 

Agsin^kya) 
Aosin{kya) y+,g+ 


(Cl) 


where x±^g^ = —(e(k) — ^) =F t 3 Cos(^^^^) and y±^g+ = (e(k) — y,) ± . We can easily check that at 

kx = 0, each block Hamiltonian owns PHS, which is the ID Kitaev model for p-wave superconductors. Based on the 
above block diagonal Hamiltonian, one can see that it can be written as 




a;+_g+(k) Aosin^kyu) 
Aosin{kya) y+,g+(k) 


( 02 ) 


in the extended Brillouin zone, where kx C [—%■, and C [—fj f-]- At momenta kx = 0,^, the Hamiltonian 
respects PHS. Thus, we can define a Z 2 topological invariant Vg^, which is expressed as (—l)'"s+ = sign{\t 2 \ — 
|^_)_|) 30i 40,42 , .^yjjere = g — uiq — ti — t^cos{^) at = 0 and /2-|. = /i — mg — G + tscos^^) at kx = When 
we fold the extended BZ, kx = 0, are all mapped at kx = 0. Therefore, possible Majorana zero modes appear at 
kx = 0. It should be pointed out that the topological invariants are different at kx = 0,^ by choosing appropriate 
parameters tg and (p. 

For the case with the G_ symmetry, G_(k) = {cos{^^)t^ 0 CTi + sin{^^)Tz 0 0 - 2 ). The eigenvalues are 

. The eigenvectors are vi^- = (0,0, , 1)^, V 2 - = (—,1,0,0)^, = (0, 0, —, 1)^ and V 2 ,+ = 


TABLE HI. Parameters for the emergence of edge modes in Fig. 2(c) and (d) for G± configurations in the main text. 



mo 


t2 

U 


<!> 

Aq 

G+ 

1.5 

- 0.5 

-3 

-1 

0 

O.Itt 

2 

G 

1.5 

- 0.5 

-3 

-0.5 

0 

O.Itt 

2 
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• kx o. rr ■ fere a 

,1,0,0)^, where corresponds to eigenvalues =Fe*^“. On the above eigenvectors, our BdG Hamiltonian can 
be expressed in a block diagonal matrix, which reads 


HBdG,g- 


X-,g- 

—Aosin^kya) 

0 

0 

Aosin{kya) 

y-,g- 

0 

0 

0 

0 

X+,g- 

—Aosin{kya) 

0 

0 

—Aosin{kya) 

y+,g- 


(C3) 


where x±,g- = —(e(k) — /r) ± t 3 Cos{ a^d y±,g- = (e(k) — y.) ± t 3 Cos{ can easily check that at 

kx = ^, each block Hamiltonian owns PHS, which is the ID Kitaev model for p-wave superconductors. 

The Hamiltonian for G_ in the extended BZ can be written as 




®+.s-(k) 

—Aosin{kya) 


~Aosin{kya) 

2/+.s-(k) 


(C4) 


where kx C and ky C At momentum kx = if j the Hamiltonian respects PHS. The Z 2 topological 

invariant Vg- is defined as (—= sign{\t 2 \ — |/r+l)j where fi+ = g — mo — H T t 3 sin{^) at momenta kx = if. 
When we fold the extended BZ, kx = i- are mapped at kx = i-. Therefore, possible Majorana zero modes appear 
at kx = i-. 

Another feature that can be extracted from the extended BZ is that the pairing of Cooper pairs is between two 
electrons with the momenta {k,-k) for the case of the G+ symmetry, making Majorana zero modes appear at kx = 0, 
while the pairing is between two electrons with the momenta (fc,Q — k) with Q = (^,0), which is known as the t] 
pairing, for the case of the G_ symmetry. For the 77 pairing case, the Majorana zero modes emerge at fca, = ^ = f j 
where kx = — + Q. We do the Fourier transformation of the gap functions from momentum space to real space for 

only X direction. For normal paring functions in the model of extended BZ, A{6rx)norm = = 

jfJ2i^{cG,kyCi^-Sr^,-ky), where A(k) = (ckc_k) and Ck = are used. This result indicates that 

the pairings along y direction are the same for A and B sites if 5rx = 0. On the other hand, for rj paring functions 


A(5rx)g — A(Q kx-,ky)e 


ikxSrx 


~ N Si 


a ^Q'i> 


{^ix,ky^ix—Srx, — ky) jy ^ {^ix,kyCix—Srx, — ky)i where 


integers n = {0,..., 2N} denote index of sites along x direction. Thus, the 77 pairings along y direction are of opposite 
signs for A and B sites if Svx = 0. 


Appendix D: Topological classification for two-dimensional superconductors with nonsymmorphic crystalline 

symmetry 


1. Superconductors without time reversal symmetry 

The BdG Hamiltonian for both spin-i and spinless superconductors without TRS belongs to symmetry class 
D. Table IIVI lists topological classifications for two-dimensional(2D) TNGSc (glide plane symmetry here) in the 
symmetry class D, in which the topological classification is Z in 2D. When the glide plane symmetry is included, the 
topological invariant can be described by two ID topological invariants at the momentum lines k • r = 0, We list 
the corresponding symmetry classes and also their topological invariants along these momentum lines for Gy in the 
last two rows of Table IlYl In the following subsections, we analyze the Hamiltonian symmetry classes at momentum 
lines k ■ r = 0, ^ to obtain topological classification^. In this section, we focus on the case with a zero Ghern number 
in the whole 2D Brillouin zone and the Z 2 topological invariants in the two glide parity subspaces would be the same, 
leading to only one Z 2 invariant at momentum lines k • r = 0, The case with non-zero Chern number in the whole 
2D Brillouin zone is discussed in Sec. V(D) of the Supplemental Material and also in Ref. IM 


a. class D: G+ for spinless superconductors 

For GG+(k)G“^ = G+(—k), we select a state ipk with glide parity with S = ±1. Its PH partner 7 /!_k has 

glide parity Thus, at the momentum line k • r = 0, PHS exists in each glide parity subspace. Therefore, 

the symmetry class in each glide parity subspace is D and topological invariant in ID is Z 2 . At k • r = ^, PHS 
interchanges states in opposite glide parity subspaces. Thus, the symmetry class in each glide parity subspace is A 
and there is no topological classification in ID. 
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TABLE IV. topological classification of 2D nonsymmorphic crystalline superconductors (glide plane symmetry here) without 
TRS. The class of Hamiltonian and its corresponding topological classification are listed in one block. NoGS in the hrst column 
represents a system has no glide plane symmetry. stands for no topological classihcation. 



spinless 

spin-^ 

NoGS(2D) 

D, Z 

D, Z 


k ■ r = 0 

k.r=? 

k ■ r = 0 

II 

b- 

G+ 

D, Z 2 

A,- 

A,- 

D, Z 2 

G- 

A,- 

D, Z 2 

D, Z 2 

A,- 


b. class D: G- for spinless superconductors 

For C'G_(k)C'“^ = —G_(—k), we consider a state ipk with glide parity with <5 = ±1. Its PH partner has 
glide parity —Thus, at k • r = 0, PHS interchanges states in two glide parity subspaces and the symmetry 
class in each glide parity subspace is A without topological classihcation in ID. At k • r = I-, there is PHS in each 
glide parity subspace, leading to a Z 2 classihcation in each parity subspace. 


c. class D: G+ for spin-^ superconductors 


For GG+(k)G“^ = G+(—k), the PH partner i/j-k of a state tpk with glide parity (<5 = ±i) has glide parity 

( 5 *e“®k r^ Thus, at k • T = 0, PHS interchanges states in two glide parity subspaces, giving rise to symmetry class A 
without topological classihcation in ID in each glide parity subspace. At k • r = ^, PHS exists in each glide parity 
subspace, leading to the symmetry class D and Z 2 classihcation in ID in each glide parity subspace. 


d. class D: G- for spin-^ superconductors 


For GG_(k)G“^ = —G_(—k), the PHS partner ■i/j-k of ifu with glide parity {S = ±i) has glide parity 

—Thus, at k-T = 0, PHS exists in each glide parity subspace, leading to the symmetry class D in each glide 
parity subspace and the corresponding Z 2 classihcation in ID. At k • r = ^, similar analysis suggests symmetry class 
A and no topological classihcation in each glide parity subspace. 


2. Superconductors with time reversal symmetry 

In this section, we consider the BdG Hamiltonian for spinless(spin-i) superconductors with TRS belonging to 
symmetry class BDI(DHI). We emphasize that tiem reversal symmetry always commutes with any space group 
symmetry in a physical system. Table |V] lists possible topological classihcations for 2D nonsymmorphic crystalline 
superconductors(glide plane symmetry here) in the symmetry classes BDI and DHL The topological classihcation of 
2D superconductors for class Dili is Z 2 while there is no topological classihcation for the class BDI in 2D. In our case, 
for G_, we only need to concern the two momentum lines k • r = 0, similar to the case of D symmetry class, while 
for G+, we hnd the chiral symmetry exists for any momentum. This makes the classihcation of the G+ case quite 
different from other cases. Below we will discuss each case separately. 


TABLE V. topological classification of 2D nonsymmorphic crystalline superconductors(glide plane symmetry here) with TRS. 
The class of Hamiltonian and its corresponding topological classification are listed in one block. NoGS in the first column 
represents a system has no glide plane symmetry. stands for no topological classihcation. 



spinless 

spin-4 

NoGS(2D) 

BDI, - 

Dill, Z 2 


k • r = 0 

b:|CN 

II 

1 - 

k • T 7 ^ 0 , ^ 

k • r = 0 

■s 

11 

tO|=l 

k - r / 0 , f 

G- 

AI,- 

D, Z 2 

A,- 

D, Z 2 

AH, - 

A,- 


k ■ r = 0 

k.r=^ 

k ■ r / 0 , ^ 

II 

0 

11 

k ■ r yf 0 , -7 

G+ 

BDI, Z©Z 

AHI, Z 

AHI, Z©Z 

AHI, Z 

Dill, Z 2 0 Z 2 

AHI, Z©Z 
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a. class BDI: G+ for spinless superconductors 


The TR operator is T 


^ with Qe = K for spinless systems. One can easily check that TG+(k)r ^ = 


G+(— k) and also GG+(k)G“^ = G+(—k). For a state ifk with glide parity {6 = ±1), its PHS partner tjj-u. and 

TRS partner '0_k have glide parities and respectively. Thus, at k • r = 0, PHS and TRS coexist in 

each glide parity subspace. Therefore, the symmetry class in each glide parity subspace is BDI with Z classification in 
ID. The topological invariant for the full BdG Hamiltonian is Z©Z at k-r = 0, since the two glide parity subspaces at 
momentum k ■ r = 0 are independent and there is no symmetry operation that can couple them. For the momentum 
k • r = G with 0 < G < tt/o, although both TRS and PHS do not exist in each glide parity subspace, the chiral 
symmetry, defined as the combination of TRS and PHS (H = TC), is preserved in each glide parity subspace. Thus, 
for any other momentum lines, the symmetry class in each glide parity subspace is AHI with a Z classification in ID. 
The topological invariant for the full BdG Hamiltonian is Z © Z at any momentum line 0 < k-r < 7 r/ 2 , since the 
two glide parity subspaces are independent. The TRS or PHS can relate a state at the momentum k • r = —G to a 
state with the same glide parity at the momentum k • r = G, which suggests that topological classifications in the 
momentum regime — 7 r /2 < k • r < 0 is directly determined by that in the momentum regime 0 < k • r < tt/2. At the 
momentum line k ■ r = 7 r/ 2 , chiral symmetry still exists (AHI class), and TRS or PHS relate opposite glide parity 
subspaces. As a result, the classification should be determined only by one integer Z. 


b. class BDI: G- for spinless superconductors 

In this case, we have rG_(k)T“^ = G_(—k) and GG_(k)G“^ = —G_(—k). For a state ‘ipk with glide parity 
{6 = ±1), its PHS partner and TRS partner tjj-u. have glide parities and respectively. 

Thus, at k • r = 0, PHS interchanges two states in opposite glide parity subspaces but TRS exists in each glide parity 
subspace. Therefore, the symmetry class in each glide parity subspace is AI with no topological classification in ID. 
At k • r = ^, the situation is exactly opposite and the corresponding symmetry class in each glide parity subspace is 
D with a Z 2 topological invariant. The topological invariant for the full BdG Hamiltonian is Z 2 at k • r = since 
the two glide parity subspaces at momentum k • r = ^ are related by TRS. For momentum lines k-r 7 ^ 0, the 
symmetry class is A, leading to no topological invariant in ID. 


c. class DHL G+ for spin-^ superconductors 


In this case, the TR operator is T = ^ ^ j with 0e = iaoS 2 K where a acts on different atomic sites(orbitals) 

and s acts on the spin space. Since 0eA'k(3)0(r^ = TG-|.(k)T“^ = G+(—k) and GG+(k)G“^ = G_|_(—k). 

Thus, for a state '0k with glide parity (S = ±i), its PHS partner 0_k and TRS partner 0_k have glide parities 

^♦g-ik-T j*g-jk-T^ respectively. Thus, at k ■ r = PHS, TRS and chiral symmetries coexist in each glide parity 
subspace, leading to symmetry class Dili with a Z 2 classification in ID. Because there is no symmetry relating these 
two glide parity subspaces, the topological invariant for the full BdG Hamiltonian is Z 2 © Z 2 at k • r = j. The 
chiral symmetry H exists at any other momenta in each glide parity subspace, yielding symmetry class AHI with a 
Z topological classification in ID. For the momentum line k-T = G (0<G <7r/2), the topological classification is 
Z © Z, which also determine the classification at the momentum line k ■ r = —G. It should be pointed out that in 
this momentum regime, the difference between class BDI and Dili lies in the fact that for Dili, TRS or PHS relate 
the AHI topological invariants in opposite glide parity subspaces at the momentum line k • t = G and k • r = — G 
while TRS or PHS relate the AHI topological invariants in the same glide parity subspaces for BDI class. At the 
momentum line k • r = 0, chiral symmetry exists (AHI class), and TRS or PHS changes glide parities. Therefore, the 
classification should be determined only by one integer Z. 


d. class Dill: G- for spin-^ superconductors 

In this case, rG_(k)r“^ = G_(—k) and GG_(k)G“^ = —G_(—k). The PHS partner '(/)_k and TRS partner 0_k 
of a state 0k have glide parities and respectively. Thus, at k • r = 0, TRS interchanges two 

glide parity subspaces but PHS exists in each glide parity subspace. The corresponding symmetry class is D with 
a Z 2 classification in ID in each glide parity subspace. The topological invariant for the full BdG Hamiltonian is 
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'Ll at k ■ r = 0, since the two glide parity subspaces at momentum k ■ t = 0 are related by TRS. At k • r = 
PHS interchanges two glide parity subspaces but TRS exists in each glide parity subspace, giving symmetry class All 
without any classification in ID. For momentum lines k • r 0, the symmetry class is A, leading to no topological 
invariant in ID. 


Appendix E: Models and topological invariants for TNCSc in different symmetry classes 

The toy model for TNCSc in the D class has been discussed in the main text. In this section, we will study models, 
as well as the related topological invariants, for TNCSc in other symmetry classes and show physical consequence of 
our classification in the Table IIVI and 


1. Class BDI and Dill in G_ configuration 

We first consider the case of G_ in the symmetry classes BDI and DHL In both case, the chiral symmetry 11 anti¬ 
commutes with the glide symmetry, {11, G_(k)}. This exactly corresponds to the situation discussed in Ref. Hi in 
which the model Hamiltonian belongs to AIII class with glide symmetry and the correspond topological classification 
is For the AIII class, Dirac type of edge inodes are unpinned^. For our case of BDI and Dill, we can also consider 
the model used in Shiozaki’s paper— by properly imposing additional TRS and PHS. The corresponding topological 
classification is the same {Zi), but the PHS and TRS yield the gapless point of edge modes pinned at the momentum 
k • T = 7r/2 for the BDI class and k • r = 0 for the Dili class. Topological invariants can also be defined in a similar 
manner as in Ref—. 


2. Class BDI in G+ configuration 


In this subsection, we will study a 2D model with G+ symmetry in BDI class. The model reads 
HbdiW =(e(k) - lijTsao + tscos ( —jracri + facosf —)s*n( —)t3 (T2 


,kxa^ 


,kxa. 


-Aosin{kya)T2ao - Aisin{-^)cos{-^)T2ai - Aism^(^^)T2fT2, 


(El) 


where e(k) = toq + tiCos{kxa) + t 2 Cos{kya), and and (7^ are Pauli matrices acting in the Nambu and sublattice 
space, respectively. The TRS and PHS operators read Tbdi = tqUqK and Cbdi = Tia^K. One can easily check that 
TbdiHbdi{^)Tbdi = Hbdi{-'^) and GsD/i/sm(k)G^{ 5 j = -Hbdi{-'^)- The chiral symmetry operator is Hbd/ = 
Tbdi x Cbdi = Dfo and IVbdiHbdi{}T)IV^\)i = —HbdiOT)- The glide plane symmetry operator G+ is the same as 


that in the main text, which reads G+,BD/(k) = ( 


with D]s_{g) = e* 2“ (cos(^)cti -I- sin{^)a2)- 


0 

V 0 D\{g) J . 2 2 

As we have demonstrated previously, the chiral symmetry exists for all momenta kx in each glide parity subspace, 
leading to BDI class at k • r = 0 and AIII class otherwise. The topological invariant for both BDI class and AIII class 
in ID is described by a winding number— 


“ 27]-j dkyTr[Q^^ (k)Vfcj^(5c± (^^)] 


(E2) 


for the subspace with glide parity = ±e 2 . The integral is applied along a closed loop L at momentum kx 
in the Brillouin zone. Due to chiral symmetry existing in one glide parity subspace, we can always find a unitary 

matrix U to transform the Hamiltonian into an off-block-diagonal for m^^i'^^ , UHBdGX±C^ ~ ( 

V ^C± ^ 

in which is the off-diagonal block. For our model, we can consider the integral loop L along the ky direction 
and regard kx as a parameter. Therefore, is a function of kx and Eq. (IE2I) can be further simplified as vc,^,kx- ~ 
§^d\ln{detQ(^^(kx,ky))], which suggests that the winding number is related to how many loops det[Qc^^{kx:ky)] 
evolves around the origin in the complex plane as ky changes from —tt to tt. Since TRS relates two states in the same 
glide subspace, one can easily prove the relation — ^C±,-fcx- 

To verify the Z topological invariant for BDI class in G_|_ configuration, we consider a semi-infinite system for the 
Hamiltonian (lEll) with an open boundary along the y direction and apply the iterative Green function method^^ to 



13 



FIG. 3. (Color online), (a) Edge DOS with zero energy flat bands appearing for BDI class in G+ configuration, (b) Illustration 

i k^ CL ik^ CL 

of winding number in glide parity = e~^ subspace, (c) Illustration of winding number in glide parity (- = —e~^ subspace. 
The red dot denotes the origin in the complex plane. 


calculate local density of states (LDOS) at the boundary, as illustrated in Fig. |3l(a), with the parameters listed in 
Table IVTl Strikingly, we find two-fold degenerate zero energy flat bands for all momenta kx from our calculations. 


TABLE VI. Parameters for G+ configuration in class BDI. 


mo 


t2 



< 

o 

<1 

BDI: G+ 1.5 

0.5 

2.5 

0.1 

0 

2 0.5 


In order to demonstrate topological property of these flat bands, we perform an analytical calculation for topological 

invariant. The BdG Hamiltonian (lEll) can be written in a block-diagonal form Hbdi = diag[H , H ■ Axo ], 

C+=e‘ 2 C-=-e 2 

where = (e(k) ± t3Cos{^^))T3 — Aosin{kya)T2 — Aisin{^^)T2. Each block has a specific glide parity ((±. By 
performing a further basis transformation, we obtain the off-diagonal block for , which is = e(k) ± 
t^cos{^^) — iAosin{kya) — iAisin{^^). The winding number for is mainly determined by Deltas and t 2 
in this case and for our choice of parameters, Ui^^(kx) = 1 for all momenta kx- We plot the evolution of det[Q(^^] 
for kxCL = 0.2, as shown in Fig. EKb) and (c), which confirm the non-trivial winding number for our model. Due to 
non-zero winding number Vi^^{kx) = 1, zero energy modes (Majorana modes) emerge for any momentum kx-, leading 
to zero energy flat bands. In addition, the relation = i^c±.-fcx due to TRS suggests that Majorana flat bands 

should be symmetric around k = 0 in each glide parity subspace. 


3. Class Dill in G+ configuration 

For NTCSc in class Dill, we consider the same lattice structure as that in the main text. The glide plane 
symmetry operator for the normal state reads D\^[g) = (cos(^|^)criS 3 -b sin{^^)cj 2 S^) on the basis Ck = 

(|H,k, t), |x4,k, j,), |i?, k, t), |S,k, j,))^ where at and Si are the Pauli matrices acting in the sublattice and spin space. 
The TRS operator is 0 = iaoS 2 K- The Hamiltonian, on such a basis Ck, is expressed as he(k) = e(k)croSo + 
tisin{kxa)a3Si — t^sin(kya)a^S2 + tzsin{ -b sw(^)cr2)( ^^y° ) + Hsin{ )(cos(^|^)gi -b 
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5 w(-^|^)g 2 )( 2 ^° ), where e(k) = mo + ticos{kxa) + t 2 Cos{kya). One can easily check that this Hamiltonian sat¬ 

isfies that {g)he{y)D\^{g) = /ie(k) and 0 “^he(k )0 = /ie(—k). 

On the basis 'l'(k) = (ck,ct3k)^ in the Nambu space, we can construct the BdG Hamiltonian in class Dili, which 
reads 


Hdiii =e(k)T3(ToSo + tisin{kxa)TQa^si - t5sin{kya)T^a3S2 

,, . i{kx-(t))a TQ+Tz k^a . k^a S 3 + sq 

+t 3 Sin{ - - -)(— - —)(cos( —)cti -k stn{—)a 2 )[ - - —) 

+t 3 Sin{ -^-)(— - —)(cos (—)cti -k Sin{—)a 2 ){ - - —) 

+t 3 Sin{ - - -)(— - —)(cos( —)cti -k stn{—)a 2 ){—-—) 

+t 3 Sin{ - - -)(— - —)(cos( —)cti -k stn{—)a 2 ){—-—) 

-|-Aosm(fcya)TiCToSi (E3) 

where r^, Ui and Si are Pauli matrices acting on the Nambu space, sublattice space and spin space, respectively. 

The glide plane symmmetry, TRS, PHS and chiral symmetry operators read G+^£)///(k) = ^ 

Tdiii = iToaoS 2 K, Com = TiaoSoK and Hdiii = TiCToS 2 (Note that we ignore the ‘i’ in the chiral symmetry 
operator, which is not essential). 

There is chiral symmetry in each glide parity subspace for any momenta, and thus, similar to the BDI case, we can 
also apply the ID winding number (El. IE2I) as the topological invariant to this model. One can find a unitary matrix 

^ with g'^(-k) = -q(k) 


(E4) 


V to transform the Hamiltonian into an off-block-diagonal fornt^L^ VHBdcy^ = ^ gt(k) 0 
and the corresponding winding number is defined as 


vc±,k^ = j^dkyTr[q^^{\ii)V 


where qQ^ is the block of the matrix q in the glide parity (± subspace. For Dili class, TRS relate states in opposite 
glide parity subspaces, and thus one can show However, for momentum line k • r = ^, we have 

^C+-kx=^/a = = -^c+,k^=^/a, indicating that = 0. This conclusion is consistent with the 

fact that the Hamiltonian in each glide parity subspace belongs to symmetry class Dili with TRS. However, due to 
the anti-symmetric nature of the matrix q, an independent Z 2 topological invariant^i can be defined as 


iEc±.fex=f = n, 


{Pf[qfjJC)]/^det[qc^{IC)]}. 


(E5) 


ikx a 

where C± denotes glide parity = ±ie^“ subspace, L is a loop with momentum ^ in the BZ, /C denotes time 
reversal invariant momenta {k^^ky) = (f ,0), (^, and Pf denotes the Pfafhan. = ±1 denotes trivial and 

nontrivial topological phases. 

The local DOS at the edge for a semi-infinite system is shown in Fig. Hla) with the parameters in Table IVHI 
Interestingly, we find two types of edge modes in the energy dispersion at the boundary. There are zero energy flat 
bands, similar to the case of BDI class, around T {kx = 0), and a Dirac type of energy dispersion with gapless point at X 
{kx = tt/o). To understand edge modes in Fig. HKa), we study topological invariants in each glide parity subspace. The 
Hamiltonian in each glide parity subspace is written as ^ = e{'k)T 3 ao±t 3 sin{ )( '^ 3 +tq ^^ (^ 0 + 0-3 ^ j. 


t3si„((^£^^)(l3^)(^o^) ± t3si„((^£+^)(l3|io)(£o^) ± i3si„(0£+^)(l3^)(^o^) + t4sin{k,a)Toai - 

t 3 sin{kya)Tza 2 + AQsin{kya)Tiai. The local DOS in each glide parity subspace are shown in Fig. UJc) and (d). 
We calculate the winding number Vk^ as a function momentum kx for Hdiiix±j S'® shown in Fig. HKb). The zero 
energy flat bands exist in the momentum regime of kx where = if and disappear in the momentum regime 

with — 0- Thus, the flat bands originate from the winding number similar to the case of class BDI. 

At the momentum where the winding number changes between 0 and I, the bulk superconducting gap closes. 

Thus, our model corresponds to a nodal superconductor. At the momentum kx = tt/o, although ii(^±Xx — ^ad 

= I, giving rise to Dirac type of edge modes. This confirms two independent topological invariants in our case. 
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FIG. 4. (Color online), (a) Edge DOS with both zero energy flat bands and helical edge modes for Dill class in G+ 

ifejc a 

conflguration. (b) Winding number in glide parity = ±ie 2 subspace as a function of k^- (c) Edge DOS in glide 

ikx a ikx o, 

parity subspace = ie 2 . (d) Edge DOS in glide parity (^_ = —ie 2 subspace. 


TABLE VII. Parameters for G+ conflguration in class DHL 


mo ti 

t2 tz 

ti 

^5 

p Ao cj) 

Dill: G+ -1 0.1 

1 1 

0.5 

0.5 

0 1 O.Itt 


4. Chern number in Class D 


Finally, we will show that the conventional Chern number can also be dehned in the whole Brillouin zone for 
nonsymmorphic superconductors in the D class, which just correspond to chiral topological superconductors. However, 
in this case, glide plane symmetry does not play a key role in the sense that Chern number is still well defined even 
if glide symmetry is broken. 

The model Hamiltonian of the normal state is the same as the model listed in the main text. The gap function 
reads A(k) = AQ{sin{kx)+isin{ky)), which belongs to G+ configuration. We perform a calculation of edge density of 
state(DOS) to illustrate the chiral edge mode explicitly by using iterative Green function methodic, as shown in Fig. 
[5l Here the open boundary is applied in the y direction and we consider the upper edge that parallel to x direction. 
The parameters we use are listed in Table IVHII From Fig. [5l we find one edge mode with positive velocity (a right 
mover), which is consistent with our calculation of Chern number n = 1. 
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FIG. 5. (Color online). Edge DOS with a chiral edge state appearing along the x direction. 


TABLE VIII. Parameters for G+ configuration in class D. 

mo t2 to fi (j) Ao 
Chern: G+ 1.5 -1 -1 -1 0 O.Itt 2 
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